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Abstract 

This  paper  presents  a  new,  innovative  procedure 
called  point  load  superposition  for  determining  the 
contact  stresses  in  mating  gear  teeth.  It  is 
believed  that  this  procedure  will  greatly  extend 
both  the  range  of  applicability  and  the  accuracy  of 
gear  contact  stress  analyses. 

Point  load  superposition  is  based  upon  funda¬ 
mental  solutions  from  the  theory  of  elasticity.  It 
is  an  iterative  numerical  procedure  which  has  dis¬ 
tinct  advantages  over  the  classical  Hertz  method, 
the  finite-element  method  (FEM),  and  over  existing 
applications  with  the  boundary  element  method 
(BEM).  Specifically,  friction  and  sliding  effects, 
which  are  either  excluded  from  or  difficult  to 
study  with  the  classical  methods,  are  routinely 
handled  with  the  new  procedure. 

The  paper  presents  the  basic  theory  and  algo¬ 
rithms.  Several  examples  are  presented.  Results 
are  consistent  with  those  of  the  classical 
theories.  Applications  with  spur  gears  are 
discussed . 

Introduction 

Gear  tooth  contact  stress  is  the  single  most 
important  factor  affecting  the  life  of  a  gear. 

Tooth  contact  stresses  also  have  a  significant 
effect  upon  tooth  deformations  and  thus  upon  gear 
performance  under  load.  However,  in  spite  of  their 
dominant  role  in  gearing  mechanics,  contact 
stresses  are  9till  not  well  understood  and  they  are 
not  accurately  quantified.  Contact  stresses  are 
not  nearly  as  well  documented  as,  for  example,  root 
and  fillet  stresses. 

One  of  the  reasons  that  gear  tooth  contact 
stresses  are  not  better  understood  is  that  the 
geometry  and  kinematics  of  meshing  gear  teeth  do 
not  satisfy  the  assumptions  required  for  the  appli¬ 
cation  of  the  classical  and  standard  Hertzion  con¬ 
tact  stress  analysis.1""3  The  Hertz  analysis 
assumes  a  contact  between  frictionless  cylinders. 

It  does  not  account  for  noncircular  geometry  such 
as  involute  geometry,  and  it  does  not  include  the 
important  friction  forces  which  occur  in  gear  con¬ 
tact  away  from  the  pitch  point. 

Standard  finite  element  procedures  are  also 
unsatisfactory  for  studying  contact  stresses. 
Contact  stresses  are  in  essence  stress  concentra¬ 
tions  which  are  locally  nonlinear.  Hence,  the 


standard  finite  element  procedures  need  to  be 
modified  through  the  use  of  gap  elements,  large 
numbers  of  elements,  and  iterative  procedures 
before  they  can  be  accurately  applied  in  the  study 
of  contact  stresses.  Such  modifications,  although 
feasible,4  are  burdensome  and  not  generally 
available . 

A  relatively  new  and  promising  procedure  for 
contact  stress  analysis  is  the  boundary  element 
method.  The  boundary  element  method  is  particu¬ 
larly  useful  in  regions  with  high  stress  grad¬ 
ients  -  as  with  contact  stresses.  However,  this 
method  is  still  under  development  and  applications 
with  contact  stresses  also  require  iterative  proce¬ 
dures  to  obtain  the  contact  area  geometry.' 

In  this  paper  we  present  a  procedure  dedicated 
strictly  to  contact  stress  analysis  with  direct 
application  in  gear  tooth  contact  stress  analysis. 
It  is  similar  to  the  boundary  element  method.  It 
does  not  require  the  restrictive  assumptions  of  the 
Hertz  analysis.  That  is,  with  this  method,  non- 
cylindrical  geometry  and  frictional  effects,  as 
encountered  in  loaded  meshing  gears,  are  readily 
accommodated.  The  method  is  based  upon  the  point¬ 
load  superposition  method  developed  by  Paul,9' 
and  previous  unpublished  work  of  Sisera  Jayasinghe 
at  the  University  of  Cincinnati. 

Historically,  contact  stresses  have  been  one 
of  the  most  important,  most  widely  studied,  and  yet 
most  difficult  problems  in  the  theory  of  elasti¬ 
city.  Contact  stress  analysis  can  be  traced  back 
to  the  work  of  H.  Hertz  in  18811'"'3  who  studied  the 
contact  deformation  of  smooth  elastic  bodies 
pressed  together  by  forces  normal  to  their  contact¬ 
ing  surfaces.  He  proposed  that  the  contact  region 
could  be  modeled  by  an  elastic  half  space  loaded 
over  an  elliptical  area. 

References  9  to  22  describe  recent  attempts  to 
use  the  finite  element  method  and  the  boundary  ele¬ 
ment  method  to  determine  contact  stresses. 

The  point-load  superposition  method  presented 
in  this  paper  has  been  outlined  in  a  paper  pre¬ 
sented  at  the  Health  Monitoring  for  Propulsion 
Systems  Conference  in  Cincinnati  (1990).‘3  The 
method  is  also  documented  in  a  recent  NASA  Con¬ 
tractor  Report  of  the  authors.24 

The  balance  of  the  paper  is  divided  into  nine 
part 3  with  the  following  part  summarizing  some  pre¬ 
liminary  results  useful  in  the  sequel.  The  next 


three  parts  present  the  basic  analysis.  Numerical 
procedures  are  then  discussed  followed  by  some 
simple  examples.  The  next  part  discusses  appli¬ 
cation  with  spur  gears.  The  final  two  parts  con¬ 
tain  a  discussion  and  concluding  remarks. 

Preliminary  Considerations 

Consider  two  elastic  bodies  in  contact.  Let 
the  contact  region  be  divided  into  a  mesh  of  ele¬ 
ments  (or  cells).  Let  the  forces  transmitted 
across  each  cell  be  represented  by  a  distribution 
of  concentrated  forces  acting  on  an  elastic  half 
space.  By  using  long-established  results  from  the 
theory  of  elasticity,  together  with  the  principle 
of  superposition,  the  stress  distribution  and  point 
displacements  due  to  a  cell  loading  are  obtained. 
Then,  by  superposing  the  results  from  each  cell  the 
resulting  stress  and  displacement  distributions 
within  the  contacting  bodies  may  be  obtained. 

To  develop  the  analysis  consider  an  elastic 
half  space  with  a  concentrated  normal  load  as  in 
Fig.  1. 

Let  P  be  a  typical  point  within  the  half 
space,  located  by  the  polar  coordinates  (r ,6)  as 
shown.  Then  the  radial  and  tangential  stresses  at 
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The  displacement  of  a  surface  point  are  then 
( 1-U2) 
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where  c  is  a  constant  determined  from  a  zero- 
displacement  reference  point. 

Next,  consider  an  elastic  half  space  with  a 
distributed  normal  and  tangent  La L  loading  over  a 
portion  of  the  surface  as  depicted  in  Fig.  2. 

Then  the  resulting  stresses  and  displacements 
at  a  typical  point  P  may  be  obtained  by  super¬ 
position  (integration)  of  the  foregoing  results. 
Following  the  analysis  of  Johnson,4"  the  stresses 
at  P  are: 
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In  rectangular  coordinates,  these  stresses 
have  the  form: 
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Using  the  stress-strain  equations  and  the 
strain-displacement  equations,  the  displacements  at 
a  point  P  on  the  surface  may  be  expressed  as:“ 


In  like  manner  the  tangential  and  normal  dis¬ 
placements  of  a  surface  point  are:'* 
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where  E  and  1/  are  the  elastic  constant  and 
Poisson's  ratio,  and  where  h  is  the  distance  to  a 
reference  point  of  zero  displacement  along  the 
2-axis. 

In  like  manner,  if  there  is  a  concentrated 
tangential  force  G,  the  stresses  at  P  are: 
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where  c^  and  c^  are  constants  to  be  determined 
from  the  displacement  of  an  arbitrary  reference 
point . 


where  6  is  now  measured  from  the  line  of  action 
of  force  G.  In  rectangular  coordinates  these 
stresses  are 
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Application:  Triangular  Shape 

Loading  Distribution 

Equations  (7)  and  (8)  may  be  used  to  determine 
the  effects  of  the  force  distribution  in  the  con¬ 
tact  region.  If  the  region  is  divided  into  a  mesh 
of  elements  (or  cells),  then  the  effects  of  the 
loadings  on  the  individual  cells  may  be  super¬ 
imposed.  Moreover,  the  general  loading  in  the  con¬ 
tact  region  may  be  represented  by  the  superimposed 
loadings  on  the  cells.  A  triangular  (that  is, 
piecewise  linear)  loading  distribution  over  a  cell 
forms  a  convenient  basis  for  modelling  the  general 
loading  distribution. 


where  and  C,  are  constants  to  be  determined 

from  the  displacements  of  a  reference  point  and 
where  A  and  B  are  defined  as: 
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The  contribution  of  a  triangular  load  on  an 
individual  cell  to  the  stress  at  a  typical  point 
P  may  be  obtained  by  letting  f(x)  and  g(x)  of 
Eqs.  (7)  and  (8)  have  the  forms: 

f  (x)  -  f  0  ( a-  |  x  |  )/a  ,  g(x)  (9) 
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Modeling  of  Contact  Force  Distribution 


To  illustrate  how  these  results  can  be 
superimposed  in  the  contact  region,  suppose  that 
the  contact  force  distribution  is  represented  by 
the  piecewise  linear  distribution,  as  in  Fig.  4. 


where  a  is  the  half-width  of  the  cell  whose 
center  is  at  the  origin,  and  where  fQ  and  go 
are  the  peak  triangular  normal  and  tangential  load¬ 
ings  as  depicted  in  Fig.  3. 

By  substituting  from  Eqs.  (9)  into  (7)  and  by 
performing  the  indicated  integrations,  the  result¬ 
ing  stresses  at  P  (Fig.  3)  are:25 


Let  the  contact  region  be  divided  into  n 
equal-width  elements  (or  cells)  as  shown.  Then, 
just  as  the  region  is  discretized  so  also  are  the 
force  and  displacement  distributions.  Let  a  be 
the  element  half-width.  Then  the  element  nodal 
coordinates  are  in  multiples  of  a.  That  is,  the 
x-coordinate  of  node  k  is: 
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where  rj(  r2,  r,  9^,  8 and  0  are  given  by 
(Fig.  3): 
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Similarly,  the  displacements  of  a  surface 
point  are: 
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Let  u  ^  and  u  _  represent  the  tangential 
and  normal  displacements  of  node  i  due  to  a  tri¬ 
angular  normal  load  centered  at  node  j.  Then  from 
Eg.  (12),  u  and  u  _  may  be  expressed  as 
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where  k  is  defined  as: 
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Equation  (15)  may  be  used  to  obtain  the  dis¬ 
placements  at  node  i  due  to  a  set  of  superposed 
triangular  lc  id  distribut  ions  centered  at  each  of 
the  nodes.  In  this  case  the  terms  c4if(  and 
djff^  are  interpreted  as  sums  over  j  .' 
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Finally,  if  there  is  a  horizontal  or 
tangential  loading,  an  analogous  analysis  leads  to 
displacements  at  node  i  as: 
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where  the  coefficients 


are  the  same 


as  those  of  Eq.  (16)  and  where  there  is  a  sum  on  j 
from  1  to  n. 

The  displacements  of  Eqs .  (15)  and  (18)  are 
measured  relating  to  a  convenient  reference  point, 
say  the  origin  of  the  coordinate  axes. 


by  the  friction  coefficient  so  that  alternatively 
sliding  occurs. 

Figure  7  depicts  the  cylinders  in  contact  with 
tangential  loading.  As  before  and  P,  are 

matching  surface  points  at  the  initiation  of  con¬ 
tact.  When  the  tangential  load  G  is  applied  the 
distant  points  Ql  and  Q2  have  rigid  body  dis¬ 
placements  d1>(  and  d2)<  relative  to  O.  If  the 
displacements  of  P;  and  P?  relative  to  0  are 
denoted  by  S  and  S,  ,  the  slip  S  between  p 

and  may  be  defined  as: 
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Contact  Analysis 


Equations  (15)  and  (18)  may  be  used  to  develop 
an  iterative  procedure  for  obtaining  a  detailed  and 
comprehensive  analysis  of  contact  stresses  and 
deformations . 

To  illustrate  this  consider  two  cylinders 
pressed  together  as  in  Fig.  5.  As  the  cylinders 
are  brought  together  they  initially  have  line  con¬ 
tact.  Then  as  they  are  pressed  tighter  the  line 
contact  develops  into  a  contact  strip. 

Figure  6  depicts  the  contact  region  after  a 
finite  deformation.  The  solid  curves  represent  the 
deformed  surface  profiles,  and  the  dashed  curves 
represent  the  undeformed,  or  original  profiles. 

Consider  two  points  Qi  and  Q,  relatively 
distant  from  the  contact  region.  Then  during  the 
compression  Q;  and  Q,  move  toward  the  contact 
region  and  toward  each  other  through  the  displace¬ 
ments  d,_  and  d  respectively. 


If  P  and  P,  are  in  the  nonslip  (stick) 
region  s  is  zero,  and  then 


where  D  is  the  global  tangential  displacement  of 
the  cylinders. 

numerical  Procedures 

Equations  (15)  and  (18),  together  with  the 
constraint  equations  of  the  foregoing  part,  may  be 
used  to  determine  the  nodal  normal  and  tangential 
forces  Ci  and  g^  A  difficulty  which  arises 
however,  is  that  the  extent  of  the  contact  region 
is  not  known  a-priori.  Hence,  an  iterative  pro¬ 
cedure  needs  to  be  developed  which  will  determine 
the  contact  region  as  well  as  the  nodal  forces. 

For  the  development  of  this  procedure,  it  is 
convenient  to  superimpose  Eqs.  (15)  and  (l£i  and 
express  them  in  the  form: 


Consider  two  typical  matching  surface  points 
P,  and  P„  separated  by  a  distance  h(x).  Then 
as  the  deformation  proceeds,  and  P,  simulta¬ 

neously  move  toward  each  other  and  deform  inward 
into  their  respective  cylinders.  Let  U  and  U,_ 
(measured  as  positive  into  the  respective  cylin¬ 
ders)  represent  the  z-axis  displacements  of  P. 
and  P.  toward  Q.  and  <2..  Then  if  P,  and  P. 
come  into  contact,  we  have 

Uu  ♦  U:z  ♦  h(X)  .  d.z  ♦  d2z  -  DZ  (19) 

where  is  the  global  compressive  displacement 

of  the  cylinders.  If  P,  and  P„  are  outside  the 
contact  region  we  have 

-  h(x)  >  D,  (20) 

Equations  (19)  and  (20)  can  be  used  to  deter¬ 
mine  the  extent  of  the  contact  region. 
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where  c  and  d  are  defined  as 
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where  c  and  d  are  influence  coeffi~**nts  of 
the  origin  (Eqs.  (16)  and  (17)).  Then  \ \  substi¬ 
tuting  from  Eq9.  (24)  and  (25)  into  Er  42Ci  we 
have 

X;C  f,  .  X;d  g  .  -b1  i,j»0 . in  (27) 


If  there  is  no  relative  horizontal  displace¬ 
ment  of  P;  and  Pjf  we  have  the  no-slip  condition 


When  the  cylinders  are  subjected  to  tangential 
loading  in  addition  to  the  normal  loading,  the 
deformation  geometry  depends  upon  the  relative  mag¬ 
nitudes  of  the  normal  and  tangential  loads.  The 
tangential  loading  causes  an  asymmetric  deformation 
pattern.  This  asymmetry,  or  distortion  is  limited 
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with  E  ,  v^.  E,,  and  t',  being  the  elastic  con¬ 
stants  of  Bodies  1  and  2. 


Similarly,  for  those  nodes  in  the  stick  region 
Eqs.  (24)  and  (25)  together  with  Eq.  (21)  lead  to 
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axial  length  and  no  tangential  loading.  The  cylin¬ 
ders  had  radii:  100  and  150  in;  elastic  modulus: 
30xl0b  Ib/in";  and  Poisson  ratio:  0.3.  The  cylin¬ 
der  surfaces  were  taken  to  be  frictionless  so  than 
a  comparison  could  be  made  with  the  classical  Hertz 
solution,  which  assumes  frictionless  contact. 

Table  1  shows  a  comparison  of  numerical 
results  and  the  Hertz  solution. 
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Finally,  for  those  nodes  in  the  slip  region  we 

have 

g  *  tfif ,  i  *  -n,  .  .  .  ,  p,  p*m, .  .  . ,  n 


(31) 


where  the  sign  is  chosen  so  that  the  force  direc¬ 
tion  ’s  opposite  to  the  direction  of  slip. 

In  Eqs.  (27),  (29),  and  (31)  there  are  2n  ♦  1 
unknown  f  and  2n+l  unknown  g^  (i  *  0,...,±n) — 

that  is  4n  ♦  2  equation.  Equations  (27),  (29), 

and  (31)  constitute  4n  +  2  equations.  However, 
the  system  is  not  complete  since  when  i  is  zero, 
Eqs.  (27)  and  (29)  are  degenerate.  Hence,  two 
additional  equations  are  needed.  These  are 
obtained  from  global  force  summation  leading  to: 

F  =  f  1Ai  and  G  -  g.Ai 


(32) 

where  Ai  is  the  element  area  associated  with  node 


The  numerical  procedure  is- then: 

(1)  Assume  a  contact  region 

(2)  Divide  the  contact  region 
into  elements  or  cells 

(3)  Assume  a  stick  region  within 
the  contact  region 

(4)  Solve  Eq9.  (27),  (29),  and 
(31)  for  the  nodal  forces 

(5)  Adjust  the  extent  of  the 
contact  region  by  deleting  nodes 
with  negative  normal  forces 

(6)  Adjust  the  stick  region  using 
Eq.  (31) 

(7)  Repeat  steps  (1)  to  (6)  until 
convergence  is  obtained 

This  procedure  is  dependent  upon  superposing 
the  forces  from  the  nodes  in  the  contact  region. 
Therefore,  the  method  is  called  point  load 
superposition. 


Next  to  examine  the  effects  of  friction  (and 
nence,  the  consequences  of  neglecting  friction),  we 
considered  two  unlubricated  cylinders  with  radii 
100  and  150  in.  made  of  steel  and  aluminum  respec¬ 
tively.  Different  materials  were  chosen  to  illus¬ 
trate  the  effects  of  material  properties  on  the 
solutions.  (The  Hertz  solution  is  based  upon  iden¬ 
tical  material  properties  of  the  contacting 
bodies.)  As  in  the  first  example,  the  cylinders 
were  pressed  together  with  a  normal  force  of  1000 
lb  per  unit  axial  length.  The  elastic  moduli  for 
steel  and  aluminum  were  taken  as:  30x10**  and 
12xl06  lb/in."  and  the  Poisson  ratios  as:  0.3  and 
0.33  respectively.  The  friction  coefficient  was 
assigned  as:  0.5  to  simulate  the  unlubricated 
surfaces.  Table  2  shows  a  comparison  of  numerical 
results  and  the  Hertz  solution.  Figure  8  depicts 
the  loading  distribution  in  the  contact  region. 

Finally,  to  examine  the  effects  of  sliding, 
the  steel/aluminum  cylinders  of  the  previous  exam¬ 
ple  were  also  loaded  tangentially  to  produce 
sliding.  Table  3  shows  a  comparison  of  the 
numerical  results  and  the  Hertzian  solution. 

Figure  9  shows  the  normal  pressure  distri¬ 
bution  obtained  from  the  Hertz  and  numerical  solu¬ 
tions.  (Observe  the  shift  in  the  numerical 
distribution  due  to  the  sliding. )  Figures  10  and 
11  show  the  shear  stress  distributions  on  the  steel 
and  aluminum  cylinders  for  the  Hertz,  friction  and 
sliding  cases. 

Application  with  Involute  Spur  Gears 

To  study  contact  stresses  in  mating  gear  teeth 
we  model  the  teeth  as  contacting  cylinders  as 
depicted  in  Figure  12.  The  radii  of  the  cylinders 
depends  upon  the  position  of  the  contact  point 
along  the  line  of  contact. 

In  Fig.  12  C(  and  C,  are  the  centers  of  the 
meshing  gears  and  P  is  the  pitch  point.  r,,  I,, 
°l/  0,,  P,  and  O  are  points  along  the  pressure 
line  which  is  tangent  to  the  base  circles  cf  the 
gears  as  shown.  I ^  and  I ,  are  at  the  inter¬ 
sections  of  the  pressure  line  and  lines  perpen¬ 
dicular  to  the  pressure  line  passing  through  C, 
and  respectively.  and  O,  are  at  the 

intersections  of  the  pressure  line  and  the  addenda 
of  the  gears.  O  is  a  typical  contact  point  along 
the  pressure  line. 

Then  the  base,  pitch,  and  addenda  circle  radii 
are  given  by 


Example  Solutions 


i c  i 1 1 1 


|C;I-| 


I  c  i  F  I 


To  illustrate  the  efficacy  of  the  method  we 
flrat  considered  two  frictionle9s,  steel  cylinders 
pressed  together  with  a  force  of  1000  lb  per  unit 


p2 


-  I  C,P| 


-  I  C.o, 


I  c:°:l 


(33) 
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where  the  subscripts  1  and  2  identify  the  associ¬ 
ated  gear,  and  the  notation  |  cx  I ” |  for  example, 

designates  the  length  of  the  line  segment  between 
C  and  II, 

Since  these  radii  are  known  for  a  given  pair 
of  meshing  gears,  it  is  convenient  to  express  the 
tooth  curvatures  (and  hence,  the  contacting  cylin¬ 
der  curvatures)  in  terms  of  these  radii.  From  the 
properties  of  the  involute  profile,  the  tooth  radii 
of  curvature  at  the  contact  point  are  simply  the 
distances  from  the  contact  point  to  I,  and  I. 
respect ively . 

Consider  the  following  contact  positions: 

1.  Contact  at  the  Pitch  Point 

In  this  case  the  radii  of  curvature  of  the 
mating  tooth  surfaces  are: 

p.  -  I  TTp  I  =  I  I  c~p  I  :  -  I  cm  I  : 

(34, 


and 

P-  *  I  ip  I  *  I  I  cTp  |  :  -  !  C7T  I  :  ]:  : 

(35, 


where  the  subscripts  1  and  2  identify  the 
associated  gear. 

2 .  Contact  on  the  Addendum  Circle  of  Gear  1 


4.  Contact  at  Other  Points 

Cases  2  and  3  are  at  the  contact  extremes  with 
Case  1  being  intermediate.  For  some  other  contact 
point,  the  radii  of  curvature  will  be  different 
than  these  cases,  but  it  can  be  calculated  using 
the  same  procedures  as  in  the  above  cases. 

The  developed  numerical  procedures  may  now 
be  used  to  calculate  the  contact  stresses  on  the 
mating  teeth.  To  illustrate  this  application,  con¬ 
sider  two  identical  spur  gears  in  mesh  with  the 
following  properties. 

Number  of  teeth  18 

Pitch  diameter  3.5433  m. 

Face  width  0.3937  in. 

Pressure  angle  20° 

Transmitted  load  808  lb 

<  tangent ial ) 

Elastic  modulus  30x10'  psi 

Poisson  ratio  0.3 


The  calculated  maximum  contact  stress  and 
maximum  shear  (Von  Mises)  stress  for  two  different 
contact  points  are  presented  in  Table  4.  The  high 
contact  stress  at  the  surface  is  primarily  hydro¬ 
static  loading  and  dees  not  normally  produce  fail¬ 
ure.  The  maximum  sheer  stress  which  occurs  beneatr. 
the  surface  can  initiate  crack  formation  which 
leads  to  pitting  fatigue  failure. 

Discussion 


In  this  case  the  radii  of  curvature  are: 

p,  -  l  =  1 1  qo:  i  2  -  i  cyr  i  i:': 

(  36  | 

p. :  *  i  I  -  l  rTT  i  -  i  ro:  i  (37> 

*  (rpl  +  rp2>  3in0  '  Pi 


The  method  employs  a  double  discretat ion :  one 
of  the  contact  region,  the  other  of  the  force 
distribution.  The  execution  of  the  method  involves 
an  iterative  procedure  to  determine  the  extent  of 
the  contact  region.  In  this  way  the  method  is 
clearly  more  laborious  than  the  classical  Hertz 
method  and  even  elementary  finite-element  meth¬ 
ods  (FEM).  In  exchange,  however,  the  method  is 
believed  to  provide  a  more  accurate  and  comprehen¬ 
sive  analysis  than  either  the  Hertz  or  FEM 
techniques . 


where  <P  is  the  pressure  angle. 

3.  Contact  on  the  Addendum  Circle  of  Gear  2 
In  this  case  the  radii  of  curvature  are: 


p.  -  I  I.o,  I  -  [  |  c:o,  I  ;  -  I  CM,  I  ;  1!': 

(38) 


and 


Pi  *  Mi1;  I  *  I  I;0;  I 
•  (tp,  *  rp;>  3in^  -  P: 


(39) 


Although  these  examples  employ  bodies  with 
simple  geometries  (cylinders),  the  method  outlined 
herein  is  not  restricted  to  simple  geometries. 
Indeed,  there  are  no  restrictions  on  the  geometry 
so  long  as  the  contacting  surfaces  are  continuous 
and  geometrically  smooth  (that  is,  with  continuous 
der ivat ives ) . 

The  examples  demonstrate  the  efficacy  of  the 
method.  They  show  that  the  numerical  results  are 
extremely  close  to  those  of  the  classical  Hertz 
solutions,  for  those  cases  where  the  Hertz  solution 
is  applicable.  Indeed,  the  stresses  are  virtually 
identical  whereas  there  is  a  slight  difference  in 
the  geometric  results  -  that  is,  the  contact  region 
width  and  the  location  of  the  maximum  shear  stress. 
(The  Hertz  solution  assumes  frictionless  static 
contact  between  identical  materials.) 
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The  second  and  third  example  demonstrate  the 
significant  effects  of  friction  and  sliding  in 
contact  analysis.  The  examples  also  show  that  the 
softer  material  has  lesser  stress. 

In  addition  it  is  seen  that  friction  and 
sliding  significantly  increase  the  shear  stresses. 
Indeed,  when  there  is  sliding  with  a  relatively 
high  friction  coefficient,  the  maximum  shear  stress 
occurs  on  the  surface  as  opposed  to  being  beneath 
the  surface  (Figs.  10  and  11). 

The  conditions  of  gear  tooth  contact  (Non- 
cylindrical  surface,  friction,  and  sliding)  do  not 
satisfy  the  assumptions  for  Hertz  theory.  The 
point  load  superposition  method  is  not  restricted 
by  these  assumptions.  The  gear  stress  results  of 
Table  4  are  believed  to  be  accurate,  limited  only 
by  the  precision  of  the  descretation. 

The  method  outlined  herein  is  developed  for  a 
two-dimensional  analysis  with  a  triangular  load  on 
the  individual  elements  or  cells.  The  triangular 
load  provides  for  a  piecewise  linear  representation 
of  the  contact  loading.  The  method  can  be  extended 
to  higher  order  load  representations  (for  example, 
the  use  of  cubic  splines)  and  to  three-dimensional 
analysis . 

conclusions 

A  point  load  superposition  procedure  based  on 
the  theory  of  elasticity  has  been  used  to  study 
contact  stresses  between  elastic  bodies,  including 
meshing  spur  gear  teeth.  This  method  is  validated 
by  comparison  with  results  from  the  classical  Hertz 
method  for  frictionless  cylinders  in  contact.  The 
following  specific  conclusions  are  reached: 

1.  A  new  method  of  contact  stress  analysis 
has  been  presented.  It  is  a  numerical  method  which 
has  fewer  restrictions  than  classical  procedures. 

2.  The  method  i3  applicable  for  a  broad  range 
of  geometries  of  the  contacting  bodies.  The  con¬ 
tacting  bodies  may  have  different  material  pro¬ 
perties.  The  contacting  surfaces  need  not  be 
smooth.  Sliding  between  the  contacting  surfaces 
is  also  permitted. 

3.  The  accuracy  of  the  method  is  demonstrated 
by  a  comparison  with  results  from  classical  methods 
for  simple  cases  where  the  classical  method  is 
appl icable . 

4.  The  established  efficacy  of  the  method 
justifies  applications  with  contacting  gear  teeth. 
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Table  1.  -  Comparison  of  Hertz  and  Numerical  Solutions  for 
Contacting  Steel  Cylinders  without  Friction 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width  of  Contact  Region 
(in.  ) 

0.136156 

0.137287 

0.83% 

Maximum  Normal  Pressure 
(psi) 

9351. 3S 

9358.90 

0.08% 

Maximum  Shear  Stress  (psi) 

2808.05 

2810.26 

0.08% 

Location  of  Maximum  Shear 
on  Z-Axis  ( in. ) 

0.05310 

0.05352 

0.79% 

Table  2.  -  Comparison  of  Hertz  Solutions  (Identical,  Frictionless 
Materials)  and  Numerical  Solutions  for  Contacting 
Steel/Aluminum  Cylinders  with  Friction  (The  numbers 
in  parenthesis  refer  to  the  aluminum  cylinder) 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width  of  Contact  Region 
(  in.  ) 

0.18896 

0.18873 

0.12% 

Maximum  Normal  Pressure 
(psi) 

6738.12 

6810.89 

1.08% 

Maximum  Shear  Stress  (psi) 

2021.43 

2081.32 

(2031.28) 

2.95% 
(0.48% ) 

Location  of  Maximum  Shear 
on  Z-Axis  ( in. ) 

0.07369 

0.07967 

(0.06461) 

8.1% 

(13%) 

Table  3.  -  Comparison  of  Hertz  Solutions  (No  Sliding)  and  Numerical  Solu- 


tions  for  Contacting/Sliding.  Steel/Aluminum  Cylinders 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width  of  Contact  Region  (in.) 

0.18896 

0.19041 

0.77% 

Maximum  Normal  Pressure  (psi) 

6738.12 

6737.87 

0.00% 

Maximum  Shear  Stress  (psi) 

2021.43 

3365.39 

66.4% 

Location  of  Maximum  Shear 
on  Z-Axis  ( in . ) 

0.07369 

0.00 

100% 

Shift  of  Center  (in.) 

0.00 

0.01123 

100% 

Table  4.  -  Gear  Tooth  Stress  Results  and  Comparisons 


Pitch  Point 

Contact 

Addendum  Contact 

Surface  Contact  Stress  (psi) 

194,480 

286,013 

Maximum  Shear  Stress  (psi) 

58,407 

85,899 

Depth  to  Maximium  Shear 

Stress  (in.) 

0.0056 

0.0038 
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Figure  1.  A  Concentrated  Normal  Force  on 
an  Elastic  Half  Space 


Figure  2.  Elastic  Half-Space  with  a  Normal  and 
Tangential  Loading  Distribution 
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Figure  3.  Triangular  Normal  and  Tangential 
Loading  on  a  Typical  Cell 
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Loa*i;ni»  Distribution 


Figure  4.  A  Piecewise  Linear  Representat ion 
of  a  Force  Distribution 


Figure  5.  Two  Elastic  Cylinders 
Being  Pressed  Together 


Figure  6. 


Contact  Region  of  Deformed  Cylinder 


Figure  7.  Contact  Region  with  Tangential  Loading 


hull  of  contact  width  from  numerical  solution 
f’„  -  maximum  normal  pressure  from  numerical  solution 

- - -  normal  pressure  distribution 

n  tangential  pressure  distribution 

(direct ions  are  indicated  hv  arrows) 


Figure  8.  Contact  Region  Loading  Distribution  for 
Contacting  Steel/Aluminum  Cylinders  With  Friction 
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«  -  hull  of  contact  width  from  llci  Iziun  contact 

I©  =•  maximum  normal  pressure  from  Hertzian  contact 
u  sliding  contact 

-  Hertzian  contact 

Figure  9*  Comparison  of  Hertz  and  Numerical  Solutions 
for  Normal  Pressure  for  Contacting/Sliding, 
Steel/Aluminum  Cylinders 
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T  maximum  shear  stress 
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Figure  10.  Shear  Stress  Distributions 
on  the  Steel  Cylinder. 
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T  maximum  shear  stress 
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Figure  11.  Shear  Stress  Distributions 
on  the  Aluminum  Cylinder 


Figure  12.  Modeling  Mating  Gear  Teeth  as  Contacting  Cylinders 
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